Populations of triple negative and hormone receptor positive HER2 negative breast tumors share immune gene profiles

In breast cancer, triple negative (TN) breast cancer has most responses to immune checkpoint inhibitor (ICI) therapy. Lymphocyte infiltrate does not impact prognosis in Hormone receptor positive HER2 negative (HR + HER2-) breast tumors and few HR + HER2− tumors respond to ICI. We contrasted immune-associated gene expression between 119 TN and 475 HR + HER2− breast tumors from The Cancer Genome Atlas (TCGA) and confirmed our findings in 299 TN and 1369 HR + HER2− breast tumors in the METABRIC database. TN and HR+ HER2− tumors grouped into immune-high or -low tumors, both subtypes were represented in the immune-high group. The largest difference between the immune-high TN and HR + HER2− tumors was TN tumors had more abundant Th1 and Th2 CD4+ T cells while HR + HER2− tumors had more abundant fibroblasts (log2FC > 0.3; p < 10×10−10). This suggests an immune-high signature is not dictated by breast cancer subtype, but fibroblast subsets associated with worse outcome were higher in the immune-high HR + HER2− tumors.


Background
The histological subtyping of breast cancer, based on hormone receptors (HR) (estrogen receptor, ER, progesterone receptor, PR) and human epidermal growth factor receptor ERBB2 (HER2) is essential to clinical management, determining prognosis, and treatment of breast cancer.[1][2][3][4] The histologic subtypes have differences in the type and amount of immune in ltrate but a spectrum of immune in ltrate exist in all breast cancer subtypes.[5] The bene t of immune checkpoint inhibitors (ICI) have been seen in metastatic cancers including renal cell, melanoma, and non-small-cell lung cancer.[5][6][7] The predictive biomarkers associated with response to ICI include PD-L1 expression, IFNγ expression, and increased somatic tumor mutations.[6,7] Also, a determinant for response to ICI therapy is the presence of an adaptive immune response within a tumor prior to treatment.[8][9][10][11] Breast cancer does not have the same level of adaptive immune in ltrate seen in melanoma and non-small cell lung cancer although the level of immune in ltrate in breast cancer is associated with clinical prognosis.[12] Breast tumors with greater than 50% tumor in ltrating lymphocytes (TIL) show improved pathologic complete response to neoadjuvant chemotherapy (p ≤ 0.01) [11] and increased intratumoral CD8 + T cells predict improved overall survival (OS) (p ≤ 0.01).[13] In the adjuvant setting, each 10% increase of TIL predicted a 15% improved disease free survival (DFS) (P-value = 0.025) and 17% improved OS (P-value = 0.023) in TN but not HR + HER2-breast cancer.[14,15] In clinical studies, response to ICIs have best been seen in metastatic TN breast cancer which had a 19% response rate to pembrolizumab as a monotherapy.[16] In contrast, metastatic HR + HER2 breast cancer either with pembrolizumab as a single agent (KEYNOTE-028) or in combination with a CDK4/6 inhibitor had an overall response rate of ~ 12%.[17,18] In neoadjuvant therapy, pembrolizumab has improved pathologic complete response (pCR) in both HR + HER2-and TN.However, while the addition of pembrolizumab improved pCR in HR + HER2-tumors from 15.6-24.3%, the addition of pembrolizumab improves pCR in TN from 51.2-64.8%.[19,20] The same magnitude of bene t of ICI therapy seen in TN breast cancer was not seen in HR + HER2-breast cancer, leading to the classi cation of HR + HER2breast cancer as "cold" tumors that do not respond to immune therapy.
Targeted exome sequencing technology has proved useful for elucidating the complex relationships between cancer and the immune system through RNA expression.For example, HTG Oncology's immune-associated gene set was used to identify ICI therapy responders across a wide variety of cancers including urothelial, colon, and lung [21,22] and has been used to identify neoadjuvant ICI response in TN breast cancer in the GeparNeuvo trial.[23] The HTG oncology panel has also been used in breast cancer clinical trials to determine changes in expression that predict treatment effect with palbociclib in PALOMA2 and PALOMA3.[24] Therefore HTG gene sets have been used to investigate both immunologic and oncologic gene expression patterns in cancer that predict clinical responses.
We questioned whether we could identify distinct immune gene expression patterns in HR + HER2-and TN breast cancers that could account for the differences in immune in ltrate between these subtypes.
Analyzing gene expression data from 594 breast cancer samples using the HTG Oncology and Immunooncology panels, 119 TN and 475 HR + HER2-from The Cancer Genome Atlas (TCGA), we found two distinct subgroups of immune-high and immune-low tumors in both subtypes.This was con rmed in 1668 HER2-tumors in the METABRIC dataset, 299 TN and 1369 HR + HER2-breast tumors.After estimating relative cellular phenotype abundances from bulk RNA sequencing, we found cancerassociated broblasts more abundant in HR + HER2-tumors, particularly higher in immune-high HR + HER2-tumors than immune-high TN tumors, while helper T cell populations were more abundant in immune-high TN tumors.These differences between immune high TN and HR + HER2-tumors may contribute to the differences in response to ICI but further work is needed to determine if CAF may provide targets to improve HR + HER2-response to ICI.

Materials and Methods
Collection of datasets from The Cancer Genome Atlas (TCGA) Normalized, bulk RNA sequencing expression data and clinical annotations were downloaded from the UCSC Xena TCGA portal for 119 TN and 475 HR + HER2-(412 ER + PR+, 63 ER + PR-, 0 ER-PR+) breast cancer patients.Therefore the TCGA analysis set included 20% TN and 80% HR + HER2-tumors.HR + HER2-and TN subtypes were de ned using the clinical guidelines from the American Society of Clinical Oncology and are publicly available from previous studies of breast tumors in TCGA.[25] Mutation data (MAF) les were downloaded from the GDC Data Portal and mutational burden was analyzed using maftools.
[26] The METABRIC gene expression and annotation was downloaded from the cBioPortal for Cancer Genomics and included 1369 HR + HER2-and 299 TN tumors.Annotations for immune subtypes (C1-C6) were collected from Thorsson et al. [27] Annotations for the broblast subtypes were collected from Bartoschek et al. [28] Oncology-and immune-associated gene sets RNA expression was reported for 20,530 genes.We divided genes into two categories: oncology-and immune-associated (Table S1A-B).2,500 were captured in the HTG Edge Seq Oncology Biomarker and 543 genes were captured in the HTG Immune assays (HTG Molecular Diagnostics, Inc., Tucson, AZ).

K-means clustering and silhouette analysis
To determine the optimal number of clusters for k-means, silhouette analysis (from the R 'cluster' and 'factoextra' packages) was performed on the onco-and immune-associated genes independently.Kmeans clustering was performed (K = 2; Hartigan and Wong algorithm, 'kmeans' from the R 'stats' package) on RNA expression values for onco-and immune-associated genes groups separately.Using Silhouette analysis identi ed k = 2 as optimal number of clusters using the HTG immune genes and oncogene sets (Fig. S3A-D).Hierarchical clustering of the most variable genes in the immune-associated gene set reveals two immune-associated gene expression patterns.
Gene set variation analysis (GSVA) and cell type frequency estimation GSVA scores were calculated to measure gene set variation relative to the variation observed for all genes.To calculate sample-wise enrichment of a gene set, we used 'GSVA' in R. [29] Cell type enrichment was estimated with the 'xCell' package in R. [30] We restricted our analysis to immune relevant cell types (Table 1).To con rm xCell enrichment scores, we calculated GSVA scores for independent gene sets for CD8 + T cells (CD8A, CD8B, CD3D) or T reg (IL2RA, CCR4, CTLA4, CD28, ICOS, FOXP3, CD5).

Differential gene expression and cell enrichment
To identify differentially expressed genes or differentially abundant cell types between two groups we performed similar tests.We used limma in R [31] to calculate the log 2 fold change (log 2 FC) which describes the differences in gene expression or cell-type frequency between the designated groups.Pvalues were corrected for multiple hypothesis testing using the "BH" (Benjamini Hochberg, controlling the false discovery rate) method.

Immune-associated gene expression is not driven by clinical subtype
To investigate RNA expression pro les of breast cancer samples, we analyzed 119 TN and 475 HR + HER2-breast cancer samples and limited our analysis to the genes included in HTG Oncology's immune and oncology gene sets.For the immune-associated genes, we began by performing k-means clustering to determine the optimal number of groups to immunologically characterize breast tumors (k = 2; see Material and Methods and Fig S3A and B).Cluster 1 had higher median log 2 expression of 543 immune genes than Cluster 2 (0.29 vs. -0.32),and thus Cluster 1 was referred to as 'immune-high' and Cluster 2 as 'immune-low' (Fig. 1A).Cluster 1 contains 276 breast tumors, 84 TN tumors and 192 HR + HER2tumors (30% of the immune-high tumors are TN and 70% of the immune-high tumors are HR + HER2-).Cluster 2 contains 318 breast tumors, 35 TN and 283 HR + HER2.When we evaluated which immuneassociated genes were differentially expressed between HR + HER2-and TN breast tumors within the immune-high cluster, GATA3 (logFC = 3.8; adj.P-value = 4x10 − 101 ) was most expressed in HR + HER2tumors (Fig. 1B).GATA3 has previously been associated with hormone receptor expression and lower tumor grade.[32] We evaluated the same TN and HR + HER2-tumors using the HTG oncology-associated genes to show these tumors genetically matched the breast cancer subtypes found by immunohistochemistry.We performed k-means clustering of 2,500 oncology-associated genes and found that, as expected, there were two major clusters: Cluster 1 is 97% HR + HER2-tumors (461 tumors) and Cluster 2 is 88% TN tumors (107 tumors) (Fig. 1C).The oncology-associated genes overexpressed more in HR + HER2-are hormone related including ESR1, ERBB4, and PGR which were more expressed in cluster 1 (log FC > 5, Pvalue < 0.001 for each) (Fig. 1D).
We also performed principal component analysis (PCA) on the oncology and immune-associated gene sets individually.For the oncology-associated genes, the principal component capturing the greatest variance (PC1: 16.4%) between tumor gene expression strongly separates TN and HR + tumors (Fig. S1A).This is illustrated by an high area under the curve (AUC) of 0.96 (Fig. S1B).However, in a PCA with immune-associated genes, we nd the rst principal component (PC1: 32%) is less associated with breast cancer subtype than the second principal component (PC2: 8%; AUCs = 0.74 vs. 0.91; Fig. S1C and  D).Thus, the histologic subtype is not the primary driver of immune-associated gene expression.When evaluating all genes available in TCGA (N = 20,530) the gene expression differences between immunehigh and low breast tumors included immune response genes including T cell signaling (e.g., IL2RG) and function (e.g., GZMB) (Fig. 2A).
To provide independent validation of the immune-high and immune-low clusters, we evaluated k-means clusters in the METABRIC dataset.The METABRIC dataset includes 1658 HER2-breast tumors: 299 TN and 1369 HR + HER2-.This dataset similarly showed two immune clusters (see Materials and Methods); an immune-high cluster 1 (240 TN and 384 HR + HER2-tumors) and an immune-low cluster 2 (59 TN and 985 HR + HER2-tumors).In this data set, 80% of TN tumors and 28% of HR + HER2-tumors were in the immune-high cluster (Fig. S4A).Like the TCGA data set, when evaluating expression of the oncologyassociated gene set, breast tumors generally clustered with their HR + HER2-or TN subtype similar to results by immunohistochemistry.There were 97% HR + HER2-and 12% TN in cluster 1 and 87% TN and 3% HR + HER2-in cluster 2 (Fig. S4B).
Using immune clusters previously de ned in 33 cancer types in TCGA [27], we evaluated if these clusters distributed both within the oncology and immune-k-means clusters (Fig. 2B).We nd the immune-high samples are enriched in C2 as compared to immune-low samples (55% vs 19%; Fig. 2B i ).Immune-high contains 21% of C1 (wound healing) where immune-low contains 44% of C1.The immune subset C5 was absent in the TCGA dataset, as seen in previous studies.[27] When comparing oncology-gene based clustering, TN tumors are enriched for IFNγ dominant (C2) tumors compared to HR + HER2-(Fig.2B ii ).This is illustrated by the difference in the proportion of C2 tumors between HR + HER2-(C2 = 31%) and TN (C2 = 55%) oncology k-means clusters.This is concordant with previous reports describing a stronger T cell presence in TN versus HR + HER2-.[12] However, we nd the cluster associated with the best overall survival (C3) is almost exclusively found in HR + HER2-tumors, with 21% of the HR + HER2tumors having a C3 immune environment compared to 2.5% of the TN tumors.Forty-six percent of the C3 tumors are immune-high and, as only six of the C3 tumors were TN, the majority of immune-high C3 tumors are HR + HER2-(N = 43, 91% of all immune-high C3 tumors).[27] Inferred cell abundance differences between TN and HR + HER2-tumors within the immune-high cluster and the immune-low cluster.
To better understand the underlying differences between TN and HR + HER2-immune-high and -low clusters, we used gene signatures for CD8 + (Fig. 3A i ) and CD4 + T reg (Fig. 3A ii ) T cells as proxy for immune cell type abundance in the tumor.As expected, GSVA scores for CD8 + T and CD4 + T reg cell signatures in the immune-high cluster were twice that of the immune-low cluster.To estimate immune cell enrichment with cell type signatures, we used xCell [30] and con rmed that xCell scores were behaving appropriately in this dataset including T cells, B cells, dendritic cells, and macrophages that were signi cantly higher in the immune-high cluster (raw P-value < 10 − 50 ; Fig. 3B).
Next, we sought to characterize the immune differences between TN and HR + HER2-within the patients of the immune-high cluster.Contrasting HR + HER2-and TN in the immune-high cluster, we found TH1 and TH2 CD4 + T cells are more abundant in TN than HR + HER2-and broblasts are more abundant in HR + HER2-tumors than TN tumors (log 2 FC > 0.3; P-value < 10 − 10 , Fig. 3C).There was no association between TH1 and TH2 in ltrate in either immune-high (ρ = 0.15; P-value = 0.1859) or immune-low (ρ = 0.07 P-value = 0.6735) TN breast cancer.However, in HR + HER2-tumors there was a positive association with increased TH1 tumors also having increased TH2 in ltrate in both immune-high (ρ = 0.29; P-value = < 0.0001) and immune-low (ρ = 0.21; P-value = 0.0005) (Fig. S5).In a similar analysis of the immune-low cluster, we observe more T h1 and T h2 CD4 + T cells in TN than in HR + HER2-tumors and more broblasts in HR + HER2-tumors than TN tumors (Fig. S2).This suggests the immune in ltrate may be impacted by the differences in broblasts between HR + HER2-and TN subtype regardless of the tumors being in the immune-high or immune-low groups.
Immune checkpoint genes are differentially expressed between TN and HR + HER2-in the immune-high cluster but not in the immune-low cluster Immune checkpoint genes LAG3, ICOS, CTLA4, PD1 (PDCD1), PDL1 (CD274), PDL2 (PDCD1LG2), and OX40 (TNFRSF4), have been shown to predict responsiveness of tumors to ICI therapy.[33,34] Testing for differential expression levels of these genes, we observed statistically signi cant higher expression in the immune-high group than the immune-low group as expected (Fig. 3D, Welch's 2 sided T-test; adj.Pvalue < 9X10 − 36 ).In the immune low group the immune checkpoint target genes LAG3, ICOS, PD-1, and PD-L1 were not statistically different between TN and HR + HER2-tumors (adjusted P-value > 0.05) but CTLA4 and OX40 were higher in TNBC than HR + HER2-(adjusted P-value p = 0.016 each).Irrespective of hormone receptor status, immune checkpoint expression was signi cantly higher in the immune high cluster than in the immune low cluster (adjusted P-value < 0.001) and in the immune high group there were no signi cant differences between HR + HER2-and TN tumors.ICI e cacy is also associated with CD8 + T cells tumor in ltration and IFN-γ gene expression patterns in melanoma [35] and lung [6].We used a gene signature that contains four genes (i.e., IFNγ, CD274, LAG3, and CXCL9; IFNγ + signature), previously used to stratify urothelial and non-small-cell lung carcinoma cancer patient response to anti-PD-L1 durvalumab treatment to evaluate the immune-high and immunelow tumors.[6]When evaluating HR + HER2-and TN breast cancer, the mean of the IFNγ + signature score was higher in TN than HR + HER2-(0.08 vs. -0.047;Welch's 2 sided T-test P-value = 5.9x10 − 13 ), but the ranges of the scores were overlapping (TN: [-0.32, 0.41]; HR+: [-0.37, 0.32]).Similarly, when evaluating the immune-high group of HR + HER2-and TN tumors, the means of the IFNγ + signature scores are different, though less statistically signi cant (Welch's 2 sided T-test, P-value = 2.2x10 − 6 ; Fig. 4A).Notably, when testing the IFNγ + signature score between immune high and low cluster tumors, immune high tumors had signi cantly higher signatures scores irrespective of the breast cancer subtype (HR + HER2-Mann-Whitney test, P-value = 6.3x10 − 62 , TN Mann-Whitney test, P-value = 1.3x10 − 18 ).
Tumor mutations generate neo-antigens and immune in ltrate can increase with an adaptive immune response to these neo-antigens increasing immune in ltrate.Tumor mutational burden also can determine response to ICI that is tumor-type agnostic.
[36] To determine whether tumor mutational load differed between breast cancer subtypes, we calculated each tumor's tumor mutational burden (TMB).When evaluating the immune high group, the TN tumors contained a signi cantly higher TMB than the HR + HER2-tumors (Mann-Whitney test, P-value = 3.2X10 − 15 ; Fig. 4B).The difference in the median mutational burden between immune-high and immune-low was also statistically signi cant (Mann-Whitney test P-value = 1.2X10 − 5 ) with the immune high tumors having more TMB than the immune low tumors.The relationship between somatic mutations and the IFNγ + signature shows a moderate correlation (Spearman's ρ = 0.33), indicating that the immune response in breast cancer may in part be driven by somatic mutation derived neoantigens.Furthermore, evaluation of other mutation types (insertion, missense, nonsense) similarly showed signi cant differences between immune-high HR + HER2-and TN and immune-low HR + HER2-and TN (all p < 0.001, data not shown).
Fibroblast population associated with worse prognosis in breast cancer are higher in HR + HER2-tumors in both the immune-high and immune-low clusters.
Many different populations of cancer-associated broblasts have been identi ed in breast cancer.Four populations that have been associated with prognosis in breast cancer in the TCGA include cell cycle related (cCAF), developmental (dCAF), extracellular matrix (mCAF), and vascular (vCAF) cancer associated broblasts.The vCAF signature was an independent prognostic indicator of developing metastatic disease.In the immune-high populations, the HR + HER2-tumors expressed higher vCAF broblast signature than the immune high TN tumors (P-value = 1.7x10 − 11 ).In the immune-low, the same trend was seen but there a less signi cant difference (P-value = 7.8x10 − 3 ) (Fig. 5A).The immune-high HR + HER2-tumors also had higher mCAF expression than the TN tumors (P-value = 2.1x10 − 13 ).This mCAF signature has also been associated with increased risk of disseminated disease (Fig. 5B).
[28] TN tumors had higher cCAF and dCAF expression than HR + HER2-tumors in both immune-high and immune-low populations (Fig. 5C and D).These broblast types are not associated with tumor outcome.

Discussion
Gene expression analysis has identi ed subtypes of human breast cancer and can be used to guide treatment [37] but the differences in immune environment between the subtypes and how to address them in breast cancer remains unknown with relatively few breast cancer patients bene tting from immune based therapies.Evaluating immune gene clustering in TN and HR + HER2-breast cancer, the breast cancers did not stratify into canonical subtypes but rather divided into immune-high and immunelow tumors containing both subtypes.These immune-high and immune-low groups could be seen in two independent datasets, the TCGA and METABRIC.[38,39] Evaluation of the immune-high TN and HR + HER2-tumors demonstrated that they did not have signi cant differences in several common signatures associated with response to ICI and that the HR + HER2-immune high tumors contained both the best and worst prognostic immune environments.The largest differences in expression between immune high TN and HR + HER2-tumors were that broblasts, speci cally mCAF and vCAF that are associated with worse prognosis in breast cancer and were highest in HR + HER2-tumors while expression of TH1 and TH2 genes were highest in TN tumors.
We compared several signatures that are associated with increased immune response and improved response to ICIs including CD8 + and T reg signatures and the IFNγ signature.[40] The IFNγ + gene signature can identify tumors that respond to anti-PD-L1 ICI durvalumab in 304 NSCLC patients and 103 urothelial cell cancers independent of PD-L1 expression.
[6] The CD8 + and Treg signatures have been associated with improved response to immune checkpoint inhibitors in non-small cell lung cancer and urothelial cancer.[6,7] However, while all of these signatures were higher in the immune-high breast tumors, there were no signi cant differences between HR + HER2-and TN immune-high tumors.In immune-low tumors both CD8+, and T reg signatures were higher in TN compared to HR + HER2-(Fig.3D).
There was no difference between CD8 + and T reg signatures between immune-high HR + HER2-and TN tumors and we found HR + HER2-and TN immune-high tumors to have similar ranges of IFNγ + signature GSVA scores (Fig. 4A).In immunogenic tumors such as melanoma or NSCLC, neo-antigens were associated with T cell mediated anti-tumor responses.Increased mutational load predicted a higher rate of response to anti-CTLA4 in melanoma and anti-PD-1 in lung cancer.[41,42] There was higher TMB in TN verses HR + HER2-tumors (P-value = 3.16X10 − 15 ).Furthermore, comparing the immune high TN and HR + HER2-tumors, there was more TMB in the TN immune high tumors (P-value = 3x10 − 5 ).This might explain the increased TH1 and TH2 immune in ltrate in the TN immune high tumors as compared to the HR + HER2-immune high tumors (Fig. 4B) explaining the increased adaptive immune reactivity in TN immune-high tumors.
Six common immune environments can be seen across 33 diverse tumor types and determine tumor prognosis and response to therapy [27], with the C3 immune environment having the best clinical outcome.In this study, there was only six immune high TN tumors with a C3 immune environment therefore the majority of immune-high C3 tumors were HR + HER2-.Most immune-high TN tumors were the C2 immune subtype (Fig. 2B i ).The C6 tumors predict the worst prognosis with increased TGF-β expression and are associated with increased tumor cell migration and activation of immune suppression from the tumor stroma in breast cancer.[43,44] Only HR + HER2-immune-high tumors had the C6 immune environment (Fig. 2B ii ).In the immune-high group, HR + HER2-breast cancer contained the best prognosis (C3) and worse prognosis (C6) tumors supporting that there are differences between immune-high TN and HER2 + HER2-tumors that may explain why they have differential prognosis and response to ICI.We observed two major differences between HR + HER2-and TN within the immune-high tumors, TN immune-high tumors were enriched for TH1 and TH2 T cells while HR + HER2-immune-high tumors were enriched for broblasts, and TN immune-high tumors had increased TMB as compared to immune-high HR + HER2-tumors (Fig. 3C, 4B).
The increased expression of CAF in HR + HER2-immune-high breast tumors may suggest a mechanism that prevents these tumors from developing an appropriate anti-tumor immune response.The interaction between the tumor and stroma is important in determining the breast cancer subtype that develops.[45,46] The tumor stroma further can determine metastatic potential [47], resistance to chemotherapy[48], and modulating the immune response to the tumor [46].When looking at breast cancer as a whole, expression of stromal and broblast components have been associated with worse prognosis including PDGFRA, PDGFRB, CXCL1, CXCL14, CD10, and CD36.[49] CAFs have been associated with mesenchymal stem cells and inducing a wound healing response associated with immunosuppression.[50] In breast cancer, CAF expression has been correlated with decreased CD8 + T cells and increased macrophages.[51] Of interest, in a recent study of human breast cancer samples, four CAF subtypes were identi ed and the vascular (vCAF) and mesenchymal (mCAF) subtypes associated with developing metastatic disease.
[28] In the TCGA dataset, both of those subtypes of CAF are signi cantly increased in the HR + HER2-breast tumors as compared to the TN breast tumors and there is higher signi cance in the immune-high population (Fig. 5A and B).In this study we were unable to evaluate whether CAFs induce immunosuppression in immune-high HR + HER2-tumors but further evaluation of CAF in patients that do not respond to ICI is needed.
There are a few limitations in this study.The dataset lacks adequate tumor stage information thus we could not investigate the relationship between our ndings and the stage of breast cancer.
Previous studies have shown that more advanced tumors have decreased immune in ltrate and increased immunosuppression.[52] Most importantly, with the lack of equivalent -omics and outcome data on ICI treated tumors, we could not evaluate the e cacy of ICIs on immune-high and -low tumors or compare ICI outcomes between HR + HER2-and TN tumors.While there have been some data showing bene t in HR + HER2-patients receiving neoadjuvant pembrolizumab with chemotherapy, the magnitude of response was 24.3% pCR compared to 64.8% pCR with addition of pembrolizumab to chemotherapy in TN patients.[19,20] Still, the biomarkers that predict response to ICI evaluated showed no signi cant differences between HR + HER2-and TN in the immune-high group: HR + HER2-and TN have similar expression of CD8 + and T reg expression while immune high TN does have higher TMB.We do nd immune-high HR + HER2-tumors have increased CAF while immune-high TN tumors have increased helper T cells and increased expression of T cell modi cation genes.Further studies will need to be performed to determine whether these differences may affect how HR + HER2-and TN immune-high tumors respond to ICI.
Transcriptomic clustering of immune-and onco-associated genes in HR+HER2-and TN breast tumors.Heatmaps show log 2 gene expression for 594 tumors (y-axis).The x-axis contains the most variably expressed 25 and 28 genes in immune-and oncology-associated gene subgroups, respectively (A,C).Heatmap annotation tracks include K-means (K=2) cluster assignments, generated from either all immune-associated or onco-associated genes.Volcano plots show relationships between log 2 fold change (FC; x-axis) and Benjamini-Hochberg adjusted P-values (y-axis) for differentially expressed genes in HR+HER2-and TN subtypes using immune-associated (B) or oncology associated (D) genes.Differentially expressed genes more highly expressed in HR+HER2-tumors are red and pink, while genes more highly expressed in TN tumors are blue and light blue.In 1B, the cutoffs are (Adj.P-value<0.05;abs(Log 2 FC)>2) for red and blue, and (Adj.P-value<0.05; 1 < abs(Log 2 FC) <2) for pink and light blue.In 1D, the cutoffs are (Adj.P-value<0.05;abs(Log 2 FC)>4) for red and blue, and (Adj.P-value<0.05; 1 < abs(Log 2 FC) <4) for pink and light blue.Supplementary Files